Backflow and dissipation during the quantum decay of a metastable Fermi liquid 
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The particle current in a metastable Fermi liquid against a first-order phase transition is cal- 
culated at zero temperature. During fluctuations of a droplet of the stable phase, in accordance 
with the conservation law, not only does an unperturbed current arise from the continuity at the 
boundary, but a backflow is induced by the density response. Quasiparticles carrying these currents 
are scattered by the boundary, yielding a dissipative backflow around the droplet. An energy of the 
hydrodynamic mass flow of the liquid and a friction force exerted on the droplet by the quasiparticles 
have been obtained in terms of a potential of their interaction with the droplet. 
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D ! I. INTRODUCTION 
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\Q . A first-order phase transition at a temperature close to zero is thought to proceed via quantum tunneling through 
an energy barrier between the metastable and stable states in the configuration space. Such a quantum decay of a 
metastable Fermi-liquid phase has attracted considerable attention since it was foundo in recent experiments of the 
phase separation from supersaturated 3 He- 4 He liquid mixtures that the observed degrees of critical supersaturation 
were roughly independent of temperatures below about 10 mK. Anotheij-|possible example is a deconfinement phase 
transition in cold, dense nuclear matter as encountered in a neutron star;0 this may take place during compression of 
| the nuclear matter, induced for instance by mass accretion onto a neutron star from the companion in a close binary 
system. As long as a homogeneous medium whose pressure is near the phase-equilibrium point adjusts adiabatically 
to fluctuations of a virtual droplet of the stable phase, the rate of such a quantum nucleation may be calculated^ by 
treating the droplet radius R as a macroscopic variable, deriving its Lagrangian from thermodynamic quantities of 
the system, and describing the tunneling motion of R in a semiclassical, reversible fashion. 

In a real system, however, such an adjustment proceeds at a finite rate. Consequently, fluctuations of a virtual 
droplet are accompanied by elementary excitations in the metastable Fermi liquid. These excitations not only renor- 
malize the mass and the potential for R, but also lead to energy dissipation. Such renormalization and dissipation 
. act to alter the tunneling probability exponentiallyJa In a hydrodynamic regime where the mean free path I of the 
excitations is much shorter than i?, the particle current is completely determined,- ly the continuity equation, and the 
energy dissipation is described in terms of the kinematic viscosity of the liquid.tj When the temperature is lowered 
O ■ enough to turn the system into a colhsionless regime (R <C I), the low- lying excitations appear as quasiparticles. They 
are scattered by the droplet, and the resulting density fluctuation gives rise to a friction force exerted on the droplet 
as well as a dissipative backflow current B Here we shall investigate how the flow pattern and the energy dissipation 
depend on coupling of the quasiparticles with the droplet. n 
Backflow around a slowly moving massive impurity in a Fermi liquid has been recently calculated by Zwerger;El 
the result has been expressed in terms of phase shifts for scattering ^off the impurity. Within linear response, the 
backflow is dipolar and proportional to the density response function.Q At long distances, however, the leading term 
is not dipolar, but is the radial contribution of at least second order in the scattering potential. It was noted that 
this radial backflow arises from the asymmetry in density around-the fixed scatterer, induced by a finite background 
currents In this paper, we extend these calculations by ZwergerEI to a spherically symmetric backflow, surrounding 
a fluctuational droplet of the stable phase, in a noninteracting, neutral Fermi gas at zero temperature. We thus find 
that the backflow current induced by the linear response has the same radial dependence as an unperturbed current, 
and that these currents are governed by the particles on the Fermi surface. The friction and the dissipative backflow 
arising from scattering of these particles off the boundary are expressed in terms of the transport cross sections. Using 
the backflow pattern and friction so expressed, we derive the kinetic energy and the rate of energy dissipation of the 
metastable gas. Generalization to the interacting case is made within the Landau theory of Fermi liquids. Finally, 
the backflow and dissipation in a charged Fermi liquid relevant to dense nuclear matter are examined. 
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II. PARTICLE CURRENT AND DISSIPATION IN NONINTERACTING SYSTEMS 



We first consider a case in which a noninteracting, three-dimensional Fermi gas, composed of a single component 
(spin 1/2, bare mass m, and electric charge 0) and having an equilibrium number density m, is in a metastable state 
against a first-order phase transition at zero temperature. In this case, a spherical droplet of the stable phase, which 
is assumed to be incompressible and hence to be of uniform number density n 2 , develops inside, under, or outside 
the potential barrier. Such development, for pressures so close to the phase-equilibrium point that the droplet size is 
macroscopic, can be reduced to fluctuations of the droplet radius R(t). Let us now assume that \R\ is far smaller than 
the velocity of sound c s in the metastable Fermi gas, i.e., the gas is nearly uniform. This also means that \R\ <C vf, 
where vf = hkp/m with kp = ("S^ni) 1 ^ is the Fermi velocity. In the absence of the interaction of the droplet 
with the medium, the droplet surface acts only as a boundary condition for the unperturbed current density j(r, t)\o- 
—n 2 Ri = j(R, t)\o — niRr. Here r and f are the length and the unit vector of the position vector r, of which the 
origin is set to be the droplet center. Under such a condition, the conservation law yields 



(m — n 2 )R I — ) r , r > R 

jM)Io= ; Vr 
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r < R 



The actual current density j(r, t) deviates from j(r, t)|o by the backflow due to the interaction of particles in the 
medium with the droplet. We describe this interaction in terms of a spherically symmetric, quasistatic potential 
V(r,t) = J dr'n cx (r' ,t)v(\r — r'Q, where n ex (r,t) = n2&[R(t) —r] is the density perturbation, and v(r) is the potential 
of interaction of a particle in the medium with that in the droplet. Within linear response in v(r), the backflow is 
induced by n ex , the Fourier transform of which is given by 



n ex (q,u>) — 2-kii 2 5(uj)4:'kR 3 



sin qR cos qR 
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in the limit R — > 0. The response is characterized by the density fluctuationjll 

n(q,w)| res = w(g)x(q,w)n ex (q,w) , 



(2) 



where v(q) is the Fourier transform of v(r), and x(q, w) is the density-response function. Let us assume that the range 
a of the interaction v is short, i.e., a <C R, and set v(q) as v(q = 0)/(l + a 2 q 2 ) for simplicity. Then, the macroscopic 
change in the number of particles near the droplet may be calculated from Eq. (2) up to 0(a/ R) as 



A^V|rc 
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where the asymptotic behavior of lim g _>o x(q, 0) = —dni/dni with the chemical potential of the gas [i\ — h k F /2m 
is used. By substituting Eq. (3) into the conservation law, the radial backflow at large distances is obtained as 
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(4) 



The proportionality of j| res to aRR reflects the fact that the change A7V| res in the number of particles due to the 
linear response behaves as aR 2 . According to whether the interaction is repulsive or attractive, the backflow current 
j| res goes in the same or the opposite direction with respect to Rr. Note that Friedel-type oscillations appearing near 
the droplet are not influential as long as R is macroscopic, i.e., kpR ^> 1. Then, the backflow of form (4) is expected 
to approximately express the flow pattern for r > R. The kinetic energy K of the metastable Fermi gas can thus be 
calculated from the current sum j mac = j|o + j|rcs as 
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It is important to note that K is directly relatecfl to the reversible motion of R, since j mac satisfies the conservation 
law onra macroscopic time scale of R/\R\. In the incompressible limit, Eq. (5) reduces to the result of Lifshitz and 
Kaganu 



2 



The macroscopic current j mac in turn causes a backflow current beyond linear response through the scattering 
of particles in the medium off the droplet. Let us assume the droplet to be inert and characterize the behavior of 
the particles by the incoming momenta k and the outgoing scattering states in V: ^fc(r) = (r|k+). By noting the 
correspondence with a local reference frame moving with the velocity v mac = j ma c/^i, the equilibrium distribution 
of the particles is obtained as a shifted Fermi distribution function, n°(ek-^. Vmac )7 with = h 2 k 2 /2m. The actual 
current density may thus be calculated up to linear order in R as 



(2tt) 3 



dk — Im 

m 



^(r)j^fc(r) 



C^fcVmac • klm 



^(r)-^(r) 



(6) 



k—kf 



where k is the unit vector of k, and dflk is the integration over the directions of k. Equation (6) shows that the flow 
is governed by the particles on the Fermi surface. At large distances, by using such an asymptotic form of ipk( r ) as 
e zk-r _|_ j fc (k • f )e lkr /r + . . . with the usual scattering amplitude fk, we obtain the backflow current <5j | curr = j j 
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where <7tr(fc) = J (iSlfc (1 — k ■ f)|/fc| 2 is the exact scattering cross section. Here we have used the asymptotic form 
e i(fcr-k r) _^ ig(fi k — fi r y2kr as well as the optical theorem. We note that j mac comes from Eq. (6) via the unperturbed 
term e lk r . We see from Eq. (7) that the backflow <5j| CU rr is at least of order V 2 and is induced by the macroscopic 
current j mac in its opposite direction. 

The physical origin of the backflow i5j| curr is analogous to that given by Landauerp for a backflow around a slowly 
moving impurity. A portion of the metastable Fermi gas, which carries the macroscopic current j mac but cannot adjust 
adiabatically to the development of the droplet, is scattered by the boundary, leading to a density fluctuation around 
the droplet. This fluctuation in turn diffuses at a velocity of vp, yielding the asymptotic backflow current written by 
Eq. (7). This current has nonvanishing divergence, since the density fluctuation is positive or negative definite. This 
is a contrast to the case of a slowly moving impurity in which the density fluctuation, being positive in front of and 
negative behind the scatterer, is distributed in such a way that the corresponding backflow has vanishing divergence. a 

We proceed to calculate, the friction force F due to momentum transfer from the medium to the boundary. This 
force may be expressed astJ 
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F originates from the density fluctuation induced by the macroscopic current j mac , since it vanishes when v mac \ r =R = 0. 
Up to O(R) and 0(a/R), we obtain 
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Here we have used the relational (k+ |<9V7<9r|k+) = (h 2 k 2 /m)er t r(fc)k, derived from the Lippmann-Schwinger equation. 
The rate of energy dissipation due to this friction may then be calculated as E = — 4irF ■ v mac | r= fl. From Eqs. (1), 
(4), and (9), we obtain 
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(10) 



d/ii ri2 — n\ V3 ' R 

with a — at I (kp)/TrR 2 . We thus find that Eq. (10) is similar to the result obtained by Burmistrov and Dubovskiii 
using a dimensional analysis based on the expression for E in the hydrodynamic regime. Their result, however, has 
left unknown the coefficient of nihkp(l — n,2 / ni) 2 R 2 R 2 , which is now expressed in terms of the potential v. 
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III. EXTENSION TO INTERACTING AND CHARGED SYSTEMS 



Short-range interaction v between particles in the medium plays a role in modifying the backflow current and 
the rate of energy dissipation calculated for a noninteracting Fermi gas. This role can be examined within the 
Landau theory of Fermi liquids ,0 which allows one to construct a quasiparticle state ip p with momentum p from the 
corresponding noninteracting state via adiabatic switching on of the interaction v. In order to obtain the distribution 
n p of the quasiparticles, it is sufficient to notice that a transport equation for the quasiparticles is identical with the 
corresponding equation for a noninteracting system, except that the velocity term is determined by the quasiparticle 
energy e p and that the force —de p /dr, exerted by other quasiparticles, emerges. The transport equation thus reads 

dn p + dnp _ dip _ drip d(e p + V) _ Q 
dt dr dp dp dr 

Let us divide n p into the equilibrium part (n p ) and the fluctuating part 5n p . For the translationally invariant system 
considered here and dn p /dt = 0, we obtain (n p ) = n°{e\ ), where e£ = h 2 k 2 /2m* with the quasiparticle 

ftP h V mac 

effective mass m* , and 5n p = (n p )(\ip p \ — 1). Note that (n p ) and Sn p , associated with the current densities j ma c and 
<y|curr, reduce to the similar expressions n°(ei n _n> ) and n°(ei ™. XlV-VfiJ 2 — 1) f° r v = 0. Consequently, jL, 

frr ft mac frl* mac f I 

j| rcs , and (5j| C urr are still described by Eqs. (1), (4), and (7) in which the scattering amplitude fk and the thermodynamic 

quantity d/ii/drii are modified by the interaction £i. Here, the interaction potential v is assumed to be unchanged 

by v. This holds if the potential v is a scalar fieldE By using the obtained current densities and distribution of the 

quasiparticles, we can calculate the kinetic energy K, the friction force F, and the rate of energy dissipation E. The 

results are given by Eqs. (5), (9), and (10) in which the scattering amplitude fk and the thermodynamic quantity 

d[i\jdri\ are those calculated for an interacting system, and the momentum of unit volume n\hhp is multiplied by 

m/m*. Note that the generalization to the interacting case mentioned above is valid for a Fermi liquid at the low 

temperatures that ensure I 3> fi£3 At higher temperatures, where quasiparticle collisions should be built into the 

transport equation as a collision integral, the current and the rate of energy dissipation remain to be examined. 

Next, with reference-to dense nuclear matter that is neutralized by a gas of electrons and muons and is metastable 

against deconnnementjj we consider a case in which a metastable Fermi liquid consists of particles of charge q±e that 

interact mainly via the short-range potential v and are embedded in a uniform neutralizing background, and a stable 

phase is immersed in the same background. Then, a droplet of the stable phase has nonvanishing charge as long as 

rii ^ 71,2. This charge is in turn screened by the charged particles in the medium. The resulting particle migration is 

characterized by the Thomas-Fermi screening length Atf = [{d^i/dnij/A-K^qie) 2 ] 1 / 2 , since \R\ <C c s , being assumed 

in the medium, is consistent with the condition Atf 3> R, appropriate for the linear Thomas-Fermi analysisEJ of 

the charge response of the medium. Let us take account of the electrostatic potential <j) of interaction between the 

particles inside and outside the droplet and the charge density perturbation p cx = q\e{u2 — ni)8[R(t) —r], in addition 

( o 

to the short-range potential v and the density perturbation n cx . Accordingly, the Coulomb- induced backflow j|res of 

linear order arises along with the backflow j| rcs given by Eq. (4). When r » Atf is satisfied, one obtains j | i- cs = — j|o 

(c) 

since perfect screening ensures cancellation between the unperturbed current j| and the backflow current j|r CS ■ At 

distances r < Atf, however, the screening action of the particles gives rise to slight inhomogeneity of order (R/\tf£ 

in the density distribution of the medium, which we have calculated within the linear Thomas-Fermi approximation.E3 

..(C) 

The resulting density distribution yields the backflow j|rcs via the continuity equation: 



J 



\W ={-! + [! -KR + 0( K 2 R 2 )}(l + K r)e- K ^}j\ Q , (12) 



with k — A T p. Up to O(kR) and O(afR), the kinetic energy K, the dissipative backflow 5j| C urrj the friction force 
F, and the dissipation rate E can be calculated by following a line of argument for a neutral Fermi liquid and by 
incorporating Eq. (12) in the macroscopic current j mac and modifications by <ft in the scattering amplitude fk- Since 
j\res ~ 0(k 2 R 2 ) at r = R, the resultant expressions for F and E are identical with Eqs. (9) and (10), respectively, 

except that riihkF is replaced by (m/m*)nihkF. On the other hand, K and <5j| curr have corrections due to the 

(c) 

expression for (5j| CU rr is given by Eq. (7) in which j mac = j|o + j|rcs + j|ros , and that for K is obtained as 



K = 2nm ni R 3 R 2 ( 1 " 2 



1 dm n 2 , .a 3 
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(13) 

In Eq. (13), the term of order kR reflects the role played by j Ires' 1 in cancelling the unperturbed current j|o at r > Atf- 
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IV. CONCLUSION 



We have evaluated the flow pattern and the friction induced during nucleation of the stable phase in a metastable 
Fermi liquid. The kinetic energy K and the dissipation rate E of the liquid, which are relevant to the calculations of 
the rate of quantum nucleation, have been expressed in terms of the interaction potential v. In order to make practical 
estimates, one needs information as to v, which is generally hard to obtain from experiments. The corresponding 
scattering potential V may have an inelastic property if the droplet undergoes excitations due to the scattering of 
quasiparticles in the medium off the boundary. It is important to recall that the results for the current and dissipation 
as obtained above are based on the assumption of the liquid being nearly incompressible. The effect of compressibility 
not only develops the nonlinear response contribution to the macroscopic current j mac , but also yields a reduction in 
the kinetic energy K via the emission of sound. The latter, which stems from the finite oj behavior of j mac , plays a 
role in increasing the nucleation rate exponentially.113 

Finally, we consider a case in which the liquid is in a superfluid state. The superfluidity has consequence in reducing 
the available quasiparticle states with momenta close to the Fermi surface and hence in weakening the friction force 
F. This effect may be estimated by changing, in the Fermi distribution in Eq. (8), the shifted quasiparticle energy 

£ k-f-v mac mto the quantity J (e£ — h 2 k 2 F /2m*) 2 + A 2 . — (m/m*)?ik-v mac + ?i 2 fc^/2m*, with the pairing gap Ak. This 
indicates that, when the macroscopic current flows slowly compared with the critical velocity that makes Ak vanish, 
i.e., |v mac | < (to* /m)Ak/hk, and the temperature is zero, the friction does not occur. This is because Cooper pairs 
sustain j mac as a supercurrent. With increasing temperature and/or flow velocity, however, quasiparticle excitations 
become prevailing, giving rise to the dissipation of energy. We thus find that the nucleation rate depends strongly on 
the temperature and the flow velocity, a feature expected to be ascertained by future helium experiments at ultralow 
temperatures. 
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